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Abstract 

We consider the incompressible Euler or Navier-Stokes (NS) equations on a 
torus T d , in the functional setting of the Sobolev spaces H™ (T rf ) of divergence 
free, zero mean vector fields on T d , for n £ (d/2+1, +oo). We present a general 
theory of approximate solutions for the Euler/NS Cauchy problem; this allows 
to infer a lower bound T c on the time of existence of the exact solution u 
analyzing a posteriori any approximate solution u a , and also to construct a 
function Jl n such that \\u(t) — u & (t)\\ n ^ 3? n (*) for all t £ [0, T c ). Both T c and 
"3i n are determined solving suitable "control inequalities", depending on the 
error of ii a ; the fully quantitative implementation of this scheme depends on 
some previous estimates of ours on the Euler/NS quadratic nonlinearity [E] 
[15] . To keep in touch with the existing literature on the subject, our results 
are compared with a setting for approximate Euler/NS solutions proposed in 
[3]. As a first application of the present framework, we consider the Galerkin 
approximate solutions of the Euler/NS Cauchy problem, with a specific initial 
datum considered in [2]: in this case our methods allow, amongst else, to prove 
global existence for the NS Cauchy problem when the viscosity is above an 
explicitly given bound. 
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1 Introduction 



In recent years, there has been some activity about approximate solutions of the 
Euler and Navier-Stokes (NS) equations, viewed as tools to infer accurate a pos- 
teriori estimates on the exact solutions. We mention, in particular: the works by 
Chernyshenko et al. [3], Dashti and Robinson [I], Robinson and Sadowski [18], and 
our papers [12] [13] [H] . The present work seats within the same research area; here 
we consider the incompressible Euler/NS equations 

du 

— = -£(u.du) + vAu + f , (1.1) 
at 

where: u = u(x,t) is the divergence free velocity field; the space variables x = 
(x s ) s=ly __ )C i belong to the torus T d (and yield the derivatives d s := d/dx s ); A := 
Ylt=i°ss is the Laplacian; (u»du) r := Ylt=i u sd s u r (r = l,...,d); £ is the Leray 
projection onto the space of divergence free vector fields; u is the viscosity coefficient, 
so that v = in the Euler case and v G (0, +oo) in the NS case; / = f(x,t) is the 
Leray projected density of external forces. The dimension d is arbitrary in the 
general setting of the paper, but we put d = 3 in a final application. 

The functional setting that we consider for Eq. (11.11) relies on the Sobolev spaces 

H" (T d ) = H" := {v:T d ^R d \ (v) = 0, divv = 0, V^A~ n v G L 2 (T d ) } , (1.2) 

with ( ) indicating the mean over T d ; for any real n, the above space is equipped 
with the inner product (v\w) n := (y/— A v\ y/— A w)l2 and with the corresponding 
norm || || n . One of the main issues in this setting is the behavior of the bilinear map 

y(v,w) := -£(vw) (1.3) 

in the above mentioned Sobolev spaces. It is well known that there are positive 
constants K n( i = K n and G n d = G n fulfilling the "basic inequality" 

p{v,w)\\ n ^K n \\v\\ n \\w\\ n+1 forne^+oo^ei^weC 1 , (1.4) 

and the so-called "Kato inequality" 

\(?(v,w)\w) n \ ^ G n \\v\\ n \\w\\ 2 n for n G (~ + 1, +oo), v G H™ , w G H^ 1 ; (1.5) 

fully quantitative upper and lower bounds on K n and G n were derived in our previous 
works pj)] [16], for reasons related to the present setting and described more precisely 
in the sequel. 

Independently of the problem to estimate K n and G n , the above two inequalities 
play a major role in the very interesting paper [3] on approximate Euler/NS solutions 
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and a posteriori estimates on exact solutions. To give an idea of the framework of 
[3] we describe a result therein, using notations closer to our setting. 

Consider the Euler/NS equation (11. ip with a specified initial condition u(x, 0) = 
Uq(x); let u a : T d x [0, T a ) — > H d be an approximate solution of this Cauchy problem. 
Given n G (d/2 + 1, +oo) (and assuming suitable regularity for u ,f, w a ), let w a 
possess the differential error estimator e n : [0,T a ) — > [0,+oo), the datum error 
estimator 5 n G [0, +oo) and the growth estimators V n ,V n+ i : [0,T a ) — > [0, +oo); 
this means that, for t G [0,T a ), 

on 

|| + £K . <9u a ) - z/An a - /) (t) ||„ < e n (t) , (1.6) 

||«a(0) - «o||n < $n , (1-7) 
IWOIU ^ ^n(t) , ||«a(0lln+l ^ ^n+l(t) (1-8) 

(with w a (t) := w a (-,t), etc.). According to [3j, Eq. (11.11) with datum u has an exact 
(strong, H™ -valued) solution u on a time interval |0,T b ) C [0, T a ), if T b (with the 
estimators for w a ) fulfills the inequality 

f n l -[ dt (G n V n (t) + K n V n+l {t)) 

5 n + / dte n (t) < — — e Jo . (1.9) 

JO (-'n-'b 

The present work aims to refine, to some extent, the approach of [3] and to apply 
it to get fully quantitative estimates on the exact solution of the Euler/NS Cauchy 
problem on T d , with some specific initial datum. Our main result can be described 
as follows: assuming suitable regularity for uo,f, w a , and intending e n ,V n ,V n+ i 
as above, suppose there is a function H n G C([0, T c ), [0, +oo)), with T c G (0,T a ], 
fulfilling the control inequalities 

^ + (G nVn + K n V n+1 )R n + GrXn + e« on [0, T c ), 0^(0) ^ 5 n (1.10) 

(with d + /dt the right upper Dini derivative, see Section |2]). Then, the solution u 
of the Euler/NS equation (11. ip with initial datum u exists (in a classical sense) on 
the time interval [0,T C ), and its distance from the approximate solution admits the 
bound 

|K*) -u a (f)||n ^ %>(t) for t G [0,T C ) . (1.11) 

Some features distinguishing our approach from [3] are the following ones. 

(i) Differently from ( 11.9j) . our control inequalities (ll.lOp depend explicitly on v and 
thus could allow a more accurate analysis of the influence of viscosity on the regu- 
larity of the Euler/NS solutions. 

(ii) Our approach promises better lower bounds on the time of existence of u. For 
example, for v > and under specific assumptions illustrated in the paper, the 
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inequalities fll.lOp have solutions Jl n with T c = +00, implying the global nature of 
the NS solution u; on the contrary, if 8 n or e n are nonzero the inequality (I1.9P cannot 
have a solution with very large T b , since the right hand side is bounded by l/(G n T b ) 
and thus vanishes for T b — > +00. 

(iii) In [3] there is not an explicit bound on the distance between u and w a , such 
as ( II. lip (however, our analysis yielding ( II. lip is greatly indebted to [3] and, in a 
sense, it mainly refines and completes a chain of inequalities for u — u & appearing 
therein) . 

(iv) The constants K n and G n in the inequalities ( II. 4p ( II. 5p are not evaluated in 
[3]. On the contrary, here we have at hand our previous results [15] [16] on these 
constants; thus, in specific applications, we can implement the control inequalities 
( II.IOP and their outcome (II. lip in a fully quantitative way. 

As an example of our approach, in the final part of the paper we consider the 
Euler/NS equations on T 3 with a specific initial datum w . Independently of the 
approach developed here, this datum has been already considered in an interesting 
paper by Behr, Necas and Wu [2J, where it is indicated as the origin of a possible 
blow-up for the Euler equations. However, in the cited work the blow-up is con- 
jectured on the grounds of a merely "experimental" analysis of a finite number of 
terms in the power series Y2t^o u i( x )t l solving formally the Euler Cauchy problem. 

In the present work, dealing with the initial datum of [2J both for v = and for 
v > 0, a different approach to the Cauchy problem is developed using the familiar 
Galerkin approximation (with a convenient set of Fourier modes), combined with our 
general setting for approximate solutions based on the control inequalities (jl.lOP ; 
in this case, the Sobolev order is n = 3, u & is the Galerkin solution and we use 
for it the required estimators, to be substituted in the control inequalities f ll.lOp 
(with the values for the constants K 3 and G3 obtained in [T5] [IS]). We search for a 
solution CR 3 fulfilling Eqs. fll.lOp as equalities (i.e., with ^ replaced by =); this gives 
rise to an ordinary Cauchy problem for which is solved very easily and reliably 
by numerical means. Admittedly, our computations are preliminary: they were 
performed using MATHEMATICA on a PC, with a fairly small set of 150 Fourier 
modes for the Galerkin approximation; we plan to develop the same approach with 
more powerful computational tools in a subsequent work. 

In a few words, our results follows: in the case v = 0, the solution 

of the control equations f ll.lOp exists on a finite time interval [0,T C ) (after which it 
blows up); so, we can grant existence for the Euler Cauchy problem on the interval 
[0,T C ), where we also have the estimate ( 11.1 ip on the H| distance between the 
exact solution u and the Galerkin approximate solution w a . (Unfortunately, T c is 
less than the blow-up time suggested in [2J for the Euler Cauchy problem, so we 
cannot disprove the conjecture of the cited paper; the situation could change using 
many more Galerkin modes, which is our aim for the future). For < v < 8, the 
situation is similar: H 3 blows up in a finite time T c , and we can grant existence for 
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the NS Cauchy problem only up to T c . On the contrary, for v > 8, our approach 
grants global existence for the NS Cauchy problem (and a bound of the type (11.1 ip 
on the full interval [0, +oo)). 

To conclude this Introduction, let us describe the organization of the paper. In 
Section [2] we present some preliminaries: these concern mainly the Sobolev spaces 
on T d , in view of their applications to the Euler/NS equations (II. ip . In Section [3] we 
define formally the Euler/NS Cauchy problem, the general notion of approximate 
solution for this problem and the related error estimators. In Section HI that con- 
tains the main theoretical results of the paper, we develop the general framework 
yielding the control inequalities ( ll.lOp . and prove the estimate (II. lip on the distance 
between the exact solution u of the Euler/NS Cauchy problem and an approximate 
solution u & (here we also give more details on the connections of the present work 
with |3j). In Section [5] we present some analytical solutions of the control inequalities 
(11.1 Op . under specific assumptions for their estimators and supposing, for simplicity, 
that the external forcing in (II. ip is zero; as anticipated, in certain cases our analyt- 
ical solutions for the control inequalities are global, thus ensuring global existence 
for the Euler/NS Cauchy problem. In Section [6] we describe the general Galerkin 
method for (11.11) ; in particular, we give error and growth estimators for the Galerkin 
approximate solutions, to be used with our control inequalities (ll.lOp . In Section 
13 we consider the Galerkin method with the initial datum of j2], both for v — 
an for v > 0. In Appendix |A] we review some comparison lemmas of the Caplygin 
type about differential inequalities; these are employed in Section @] in relation to 
the control inequalities. In Appendix (Bj for completeness we report the proof of an 
essentially known statement on the Galerkin approximants for the Euler/NS Cauchy 
problem. 

2 Preliminaries 

Dini derivatives. Consider a function 



(with T G (0, +oo)). The right, lower and upper Dini derivatives of f at any point 
to £ [0,T) are, respectively, 



f : [0, T) — > R 



1 1 — y f (t) 



(2.1) 



d+f 



dt 



(to) := liminf 

/i-)-0+ 



f(t + h) -f(t ) 
h 




oo, +oo] ; 



(2.2) 



d+f 



(to) := limsup 



f(t + h) -f(tp) 
h 




oo, +oo] . 



(2.3) 



dt 
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Of course, 

-£-(*■> « -jg-W ; (2-4) 

furthermore, the opposite function -f :i G [0,T) — f(t) is such that 

^fo) = -£(*). (,5) 

The left, lower and upper Dini derivatives — , — are defined similarly, with h — > 0~; 

dt dt 

however, left derivatives are not used in this paper. Of course, all Dini derivatives 
coincide with the usual derivative if this exists. 

Sobolev spaces of vector fields on the torus; Laplacian, Leray projection, 

and so on. We work in any space dimension d G {2,3,...} (using r, s as indices 
in {1, d}). For a, b in C d we put a»b := Ylt=i a r^r and write a for the complex 
conjugate d-tuple (a^); C d carries the inner product (a, b) h» a»b and the norm 
\a\ :— y/a»a. We often restrict the previous operations to R d . 

We consider the torus T d , i.e., the product of d copies of T := R/ (2nZ); a point of 
T d is generically written as x = (x s )s=i,...,<z- In the sequel we often refer to the space 

D(T d ) = D (2.6) 

of the real distributions on T d , and to the space 

W(T d ) = W := {v = (v r ) r= i,..., d | f r G -D for all r} . (2.7) 

Elements of D' can be interpreted as "generalized functions T d — > R d " ; in the sequel, 
we call them (distributional) vector fields on T d . D and W will be equipped with 
their weak topologies. For more details on distributions (and on the function spaces 
mentioned in the sequel) we refer, e.g., to [H]. 

Using distributional derivatives, we can give a meaning to several differential 
operators acting on vector fields, e.g. the Laplacian A : W — > W and the divergence 
div : D' — > D. Any real distribution w G D has a mean (w) := (w, l/(2ir) d ) G R 
(the right hand side in this definition indicates the action of w on the constant test 
function l/(2ir) d ). Passing to a vector field v G W, we can define componentwise 
the mean (v) G R d . 

Each vector field v G B' has a unique (weakly convergent) Fourier series expan- 
sion 

v = v ^ > v >< e Qd > e ^ x ) ■= elk ' x for ieT 'i ( 2 - 8 ) 

of course, the r-th component of Vk is (v r , e_fc) (i.e., it equals the action of v r on the 
test function e_^) . Due to the reality of v, the Fourier coefficients have the property 
Vk = V-k for all k G Z d ; one has (v) = v /(2ir) d ^ 2 . 
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In the sequel we often refer to the space of zero mean vector fields, of the diver- 
gence free (or solenoidal) vector fields and to their intersection; these are, respec- 
tively, 

D' := {v G O' | (v) = 0}, D' E := {v G W | divv = 0}, ©' E0 := D' s fl W . (2.9) 

Elements of W and D' E are characterized, respectively, by the conditions v o = and 
/c.Wfc = for all k; handling with the Fourier components of vector fields in D' , it is 
convenient to put 

Z^:=Z rf \{0}. (2.10) 

Of course, for each v G D' one has (At>) fc = — \k\ 2 Vk] this suggests to define, for any 
n G R, 

: ->■ W , v ^ y^K 1 such that ( v /Z A n w) fc := |£:| n ^ for k G Zg. (2.11) 

We denote with L 2 (T d ) the space of square integrable vector fields v : T d — > R d ; 
this is a real Hilbert space with the inner product (v\w)^ := f Td dxv(x)»w(x) 
= ^2 keZ dVk»Wk and the norm \\v\\l2 := \J (v\v) £2. For any n G R, let us consider 
the Sobolev space 

H"(T d ) = H™ := G ©' | v 73 ^ G L 2 (T d ) } = (2.12) 

= { v e W | t; = 0, J2 \ k \ 2n \ v k\ 2 < +oo } ; 
fcezjj 

this is a real Hilbert space with the inner product and the norm 

(v\w) n := (V^A n v\V^A n w) L2 = \k\ 2n v- k .w k , (2.13) 



\n-=VW)~n= \ k \ 2n \ V k\ 2 ■ (2-14) 



Clearly, n ^ n! implies H™' C H™ and || || n ^ || || n /. In this paper, we mainly fix the 
attention on the divergence free Sobolev space 

H^ (T d ) = M n so := W s n H£ = {v G K I fc-Ufc = } (2.15) 

(n G R); this is a closed subspace of H™, and thus a real Hilbert space with the 
restriction of ( | ) n . 

For m G {0, 1, 2, +00} let us consider the space C m (T d , R d ) = C m of vector 
fields v : T d — > R d of class C m in the ordinary sense; if m < +00, equip this with 
the norm 

\\v ||c m := max max \d Sl ,„ Se u r (x)\ . (2-16) 

e = 0, m, :r£ET rf 
r, si , S£ — 1, d 



There is a well known Sobolev imbedding 

C C C"\ || \\ C m ^ S mn \\ \\ n for m G {0, 1 , 2, . . . } , 71 G R, 77 > ^ + 777 , (2.17) 

involving suitable positive constants S mn d = S mn . 
For arbitrary n G R, we have 

AH™+ 2 = , Ae^+ 2 = H™ , (2.18) 

and A is continuous between the above spaces; furthermore, 

(Av\v) n = -Hn+i < -Iklln for each v G H" +2 (2.19) 

(all the above statements are made evident by the Fourier representations). 
By definition, the Leray projection is the map 

£ : W ->■ D' E , 7j £7; such that (£t;) fc = £ fc t; fc for k e Z d ; (2.20) 

here £& is the orthogonal projection of C d onto k 1 - = {a G C d | fc»a = 0} (given 
explicitly by £, k c = c — (A>c)/c/|/c| 2 and £oc = c, for k G Zq and c G C d ). For each 
real n, we have a continuous linear map 

£ \ H™ : ->■ H£ (2.21) 

(which, in fact, is the orthogonal projection of (H™, ( | ) n ) onto the subspace H" ). 

The fundamental bilinear map in the Euler/NS equations. In this subsection 
we assume n G (d/2, +oo). With this condition, we have a continuous bilinear map 
H™ x H™^" 1 — > H™, (v,w) h-> vdw, where v»dw is the vector field on T d with 
components (v»dw) r := J2^ =1 v s d s w r . By composition with the Leray projection £, 
we obtain a continuous bilinear map 

7 : H£ x H^ 1 ->■ H™ , (u, w) ^ CP(tj, w ) := -£(u.du;) . (2.22) 

This appears in the Euler/NS equations, and is referred in the sequel as the funda- 
mental bilinear map for such equations. As is known, for all v G H™ and w G H™^ 1 
of Fourier components v k and w k , CP(t>, w) has Fourier components 7(v, w) = and 

9(v, w) k = - 1 ^ v h»(k - h) £ k Wk-h (2.23) 
hez$ 

for all k G Zq, where 2 k is the already mentioned projection of C d onto k^. 

The continuity of CP is equivalent to the existence of a constant K nd = G 
(0, +oo) such that 

P{v,w)\\n < ^nlMUMU+l for U G H" , U> G H^ 1 ; (2.24) 
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we refer to this as the basic inequality about CP. With the stronger assumption 
n G (d/2 + 1, +00), it is known that there is a constant G n d = G n G (0, +00) such 
that 

\(?(v,w)\w) n \ < G n ||i;|| n ||Hln for v G EI™ , w G H£ +1 ; (2.25) 

we call this the Kato inequality, since it originates from Kato's seminal paper [BJ 
(for completeness, we mention that ('P(v,w)\w) n = —(v»dw\w) n for v, w as above). 
In our previous works [15] [IB] , we derived upper and lower bounds for the sharp 
constants in the above inequalities; throughout this paper, K n and G n are any two 
constants fulfilling Eqs. (1523]) fl2T25|) . 



3 The Cauchy problem for the Euler/NS equa- 
tions: exact and approximate solutions 

From here to the end of Section [BJ we fix any space dimension d G {2,3,...}, we 
consider the Sobolev spaces of vector fields on T d , and we choose a real number n 
such that 

n6(^ + l, +00). (3.1) 

Euler and NS equations: the Cauchy problem. Let us choose a "viscosity 
coefficient" 

z/G[0,+oo), (3.2) 

a "forcing" 

/GC([0,+oo),H^ ) (3.3) 

and an initial datum 

u G H^+ 2 . (3.4) 



3.1 Definition. The Cauchy problem for the (incompressible) fluid with viscosity 
v, initial datum Uq and forcing f is the following: 

Find u G C([0, T), H™^ 2 ) R C 1 ([0, T), H™ ) such that (3.5) 
du 

— = vAu + 9(u,u) + / , u(0) = u 
dt 

(with T G (0, +00] , depending on u). As usually, we speak of the "Euler Cauchy 
problem" if v = 0, and of the "NS Cauchy problem" if v > 0. 
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3.2 Remark, (i) The map u i— > u) sends continuously H™^ 1 into H™ , while A 
sends continuously HI™+ 2 into H™ ; this explains the appearing of HI™+ 2 in the previous 
definition, at least in the NS case. In the Euler case v = 0, A is absent from (13. 5ft ; 
so, in the previous definition and in the subsequent theoretical developments one 
could systematically replace H™+ 2 with H"^ 1 ; this is not done just to avoid tedious 
distinctions between the Euler and the NS case. 

(ii) Consider any function u G C([0, T), H" +2 ). By the Sobolev imbedding (I2.17p . 
one has 

u G C([0, T), C m ) for any m G {0, 1, 2, ...} such that n + 2 > ^ + m 

(in particular, for m = 0, 1, 2, 3) . (3.6) 



The following results are well known, and reported for completeness. 

3.3 Proposition. For any z/, / and u as above, (i)-(iii) hold. 

(i) The Cauchy problem A3.5\) has a unique maximal (i.e., nonextendable) solution, 
hereafter denote with u, with a suitable domain [0, T) (0 < T ^ +00/ All the other 
solutions of ([Op are restrictions of u. 

(ii) u has the property ( E 6]) . Furthermore, if v > and the forcing (x,t) h-> 
f(t)(x) = f\x,t) is C°° from T d x (0, T) to R d , the function (x,t) ^ u(t)(x) = 
u(x,t) is C°° as well from T d x (0,T) to R d . 

(Hi) For any t G [0, T) consider the vorticity matrix oo(t), of elements u rs (t) := 
d r u s {t) — d s u r (t) (r,s = l,...,d); let ||u;(£)||co := max rjS= i j ... j dmax T d |o; rs (t)|. // 
T < +oo ; one has 

limsup 1 1 07(^)11(70 = +00 ; (3.7) 

this implies 

limsup ||M(t)|| n = +00 . (3.8) 
t->T- 

Proof. (i) See [7]. 

(ii) The property (13.61) holds because u G C([0, T), H"+ 2 ). For v > and / of class 
C°° on T d x (0,T), the same smoothness property is granted for u by Theorem 6.1 
of [5] (with the domain D considered therein replaced by T d ). 

(iii) Let T < +00. Eq. (13.71) is the celebrated Beale-Kato-Majda blow-up criterion 
PP (see also [S]). To prove (13. 8p we note that, for each t G [0,T), 

||co>(t)||co ^ const. ||w(£) ll^i ^ const. ||u(t)|| n ; (3.9) 

here the first inequality follows from the definition of u(t) in terms of the derivatives 
d r u s (t), and the second one from the Sobolev imbedding (I2.17P (with m = 1). The 
relations (13. 7p and (13.91) immediately give Eq. (13. 8p . □ 
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Approximate solutions. Our treatment uses systematically the present terminol- 
ogy 

3.4 Definition. An approximate solution of the problem A3. 5\) is any map u a G 

C([0,T a ),H™+ 2 ) n C^fOjTa),!^) (with T a G (0,+oo]J. Given such a function, we 
stipulate (i) (ii). 

(i) The differential error ofu a is 

du 

e(u a ):=^-uAu a -?(u a ,u a )-f G C([0, T a ), U n J ; (3.10) 
the datum error is 

u a (0) - M G H^| 2 . (3.11) 
^zzj Lei m G R, m ^ n . ^4 differential error estimator of order m for u & is a function 

e m G C([0,T a ), [0,+oo)) such that ||e(u a )(t)IL ^ e m (t) forte [0, T a ) . (3.12) 
Let m G R ; m ^ n + 2. ^4 datum error estimator of order m for w a a rea/ number 
5 m G [0, +oo) swc/i that ||w a (0) — Uo\\ m ^ 8 m ; (3.13) 
a growth estimator of order m for u a is a function 

V m G C([0,T a ), [0,+oo)) snc/i tfja* \\u a (t)\\ m ^ V rn (t) forte [0,T a ) . (3.14) 

In particular the function e m (t) := ||e(u a )(i)|| m , the number 5 m := ||« a (0) — u o\\m and 
the function V m (t) := ||it a (t)|| m will be called the tautological estimators of order m 
for the differential error, the datum error and the growth of u & . 

4 Main theorems about approximate solutions 

Assumptions and notations. Throughout this section we fix a viscosity coeffi- 
cient, a forcing and an initial datum as in Eqs. ( I3.2p -( 13~4"|) (recalling the condition 
(13. ip n > d/2 + 1). We consider for the Cauchy problem (13.51) an approximate 
solution 

Ha GC([fl ) T a ) ) C 2 )nC 1 ([fl,T a ) ) C) ; (4.1) 

in the sequel e n G C([0, T a ), [0, +oo) and 5 n G [0, +oo) are differential error and 
datum error estimators of order n for w a , while V n ,V n+1 G C([0,T a ), [0, +oo)) are 
growth estimators of orders n, n+ 1 (see Definition 13 ,4p . Finally, we denote with 

ueC([0,T),M^)nC\[0,T),m n J (4.2) 

the maximal solution of (13. 5p (typically unknown, as well as T). 
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Some lemmas. For the sake of brevity, we put 

T w :=min(T a ,T) , w := u - u a e C([0, T w ), H«+ 2 ) n C\[0, T w ), H« ) . (4.3) 

Furthermore, we introduce the function 

W n : [0, r tt ) [0, +oo) , W n (t) := ||w(t) ||„ ; (4.4) 

this is clearly continuous, and C 1 in a neighborhood of any instant to £ [0,T W ) 
such that w(to) 7^ 0. In the sequel we often consider the right, upper Dini derivative 
d + W n /dt (see Eq. (12. 3p ). that is just the ordinary derivative at any to with w(to) ^ 0. 

The forthcoming two lemmas review and partially refine some results of [3] (after 
adaptation to our slightly different setting; for example, Dini derivatives are not even 
mentioned in [3]). 

4.1 Lemma. One has 

^ = vAw + ?(u a ,w) + ?(w,u a ) + ?(w,w) - e(u a ) on[0,T w ), (4.5) 

iu(0) = u - u a (0) . (4.6) 
Proof. By definition of the differential error e(w a ), we have 

du 

— ^ = vAu a + y(« a , u a ) + / + e(«a) ; (4.7) 

of course, du/dt = vAu + CP(m, u) + f whence, writing w = w a + w, 

du a dw A . , ^. . 

~dt + ltt = ^ a + w ) + T ( Ma + ^ M a + ^) + / 

= vAu a + uAw + ?(u a , u a ) + 7(u a , w) + ?(w, u a ) + 7(w, w) + f . (4.8) 
Subtracting Eq. (14. 7p from ( 14. 8 p we obtain Eq. ( 14. 5ft ; Eq. (14. 6 p is obvious. □ 



4.2 Lemma. Consider the above mentioned estimators e n ,5 n ,V n ,V n+ i, and the 
function W„ defined in Eq. ^4-4ty - Then 

< -isW n + (G n V n + K n V n+l )W n + G n W 2 n + e n everywhere on [0,T W ), (4.9) 

W n (0) < 5 n . (4.10) 
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Proof. We proceed in three steps. 

Step 1. Verification of Eq. ( 14.9)) at a time to such that w(to) ^ 0. In this case, 
w(t) ^ for all t in some interval I 3 to', W n is nonzero and C 1 on /. In the same 
interval, we have: 

d + W n dW n 1 dW 2 n 1 d, . . 1 ,dw. . , A , 

dt dt 2 W n at 2 W n dt w n dt 

^Tr( u (Aw\w) n + (?{u a , w)\w) n + (?(w , u a )\w) n + w)|w) n - (e(u a )|u>), 
£□3 W n V 

On the other hand, we have the following inequalities: 

(Aw\w) n ^-\\w\\l = -W 2 n ; (4.12) 

(y(« a ,w)H„< |(?(w a ,w)H n | ^ cgKiWMin < gaK ; ( 4 - 13 ) 

C?{w,u a )\w) n ^ \C?{w,u a )\w) n \ ^ p(w,u*)\\ n \\w\\ n (4.14) 
< K n \\u a \\ n+l \\w\\ 2 n ^ K n V n+l W 2 n ; 

(?(w,w)\w) n ^\(?(w,w)\w) n \ ^ G n \\w\\ 3 n = G n W 3 n ; (4.15) 

- (e(u a )\w) n ^ |(e(« a )|w) ri | < ||e(it a )IUMIn ^ e n W n ; (4.16) 

inserting the relations (14. 12J1 — (14.16}) into Eq. ( 14.111) . we obtain the inequality ( 14.9ft 
at all times t & I and, in particular, at time to- 

Step 2. Verihcation of Eq. ( 14.9)) at a time to such that w(to) = 0. At this instant 
we have W n (t ) = 0, and the relation (14.91) to be proved becomes 

^p(to) ^ e n (t ) . (4.17) 

To go on we recall that, given a C 1 function z : [0, T z ) — > E with values in a Banach 
space E equipped with the norm || ||, we have <i + ||z||/<it ^ ||dz/<it|| at all times, 
including instants when z vanishes (see, e.g., [T7] and references therein); if this 
result is applied to the function w with norm \\w\\ n = W n , we get 

d + VJ n / \ ,, dw , ... , . , 

— W < |l^(io)|l„ . (4.18) 

Now, we evaluate (dw/dt)(t ) via Eq. (14.5)) . taking into account that u>(t ) = 0; this 
gives 

||^(t )||„ = ||e(« a )(to)||n ^ e n (t ) , (4.19) 
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and Eqs. ffl~Tgj) ( KWf yield the thesis f HTTj) . 

Step 3. Verification of Eq. H.10\) . In fact, the definitions of W n , w and 5 n give 

W n (0) = K0)-« a (0)|| . (4.20) 
Now, the proof is concluded. □ 

To go on, we need a comparison lemma of the Caplygin type. This is a variant of 
some known results [H] [11] : for more details, see Appendix 1X1 

4.3 Lemma. Let 7 G (0, +oo]; consider a function ip G C(Rx [0, T), R), (s, £) >-)■ 
■0(s, t) possessing partial derivative d^/ds G C(R x [0, T), R); furthermore, let so G 
R. Suppose that W, ft G C([0, T), R) are functions such that 

f + -(t) < ^(W(t),t) /or all t G [0,T), W(0) < s , (4.21) 



(ft 



(t) ^ V *) for all t G [0, T), 31(0) ^ s . (4.22) 



dt 
TTien 

W(t) < 0l{t) for all t G [0, T) . (4.23) 



Hereafter, interesting conclusions will arise from the combination of the above com- 
parison result with Lemma 14.21 

The control inequalities and the main theorem. We begin with a definition, 
which is followed by the main result of the section. 

4.4 Definition. Consider a function % n G C([0,T c ), [0, +oo)) ; with T c G (0, T a ]. 
This function is said to fulfill the control inequalities (with respect to the estimators 
t-mb'm'Dm'Dn+l) if 

^ -v% t + {G n V n + K n V n+1 )R n + GrXn + e n everywhere on [0, T c ), (4.24) 

ttn(0) ^ 5 n . (4.25) 

4.5 Proposition. Suppose there is a function $l n G C([0, T c ), [0, +oo)) fulfilling 
the control inequalities, and consider the maximal solution u of the Euler/NS Cauchy 
problem $3.5\) . Then, the existence time T of u is such that 

T^T C ; (4.26) 

furthermore, 

\\u(t) - u a (t)\\ n <: R n (t) forte [0,T C ) . (4.27) 
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Proof. Let us employ the following notation, already used in the previous subsec- 
tions: T w : = min(T a ,T), w := u — w a and W n : t h- > ||iu(t)|| n (see Eqs. (14.31) (14.41) ; 
w and W„ have domain [0, T w )). We further define T := min(T c ,T) (and note that 
T c ^ T a implies T ^ T w ). For the moment we have not yet proved (I4.26p . so we do 
not know whether 7 equals T c , or not. We proceed in three steps. 
Step 1. One has 

\\u(t) -Ua(t)||n < for t G [0,7) . (4.28) 

Consider the functions W n and "Jl n on [0, 7). Due to Lemma [4.21 and to the control 
inequalities, these fulfill the relations 

*±£±(t) < ^TT^t) < ip(Wn(t),t) for all t G [0,7), W n (0) ^ 5„, ; 
at at 

d +< R 

" (t) ^V(3?n(*),t) for allt G [0,T), ft n (0) , 



where ^ : [0, T) x R -> R, ^(s, t) := -us + {G n V n {t) + K n V n+l (t))s + G n s 2 + e„(t). 
Therefore, by Lemma H~3l we have W n (t) ^ %i(t) for t G [0,T); this is just the 
relation (1428]) . 

Step 2. One has the relation (14. 26]) T ^ T c . To prove this, we assume T < T c 
and try to infer a contradiction. To this purpose we note that, according to (13. 8ft . 

limsup ||w(t)|| n = +oo . 

On the other hand, recalling Eq. (I4.28|) we can write ||«(t)|| n ^ || M a(t)|| n +||w(t) — 
Ua.(t)\\n ^ ||^a(t)||n + ^n(t) fo r t G [0,T); from here and from the continuity of 
||« a ( )||„, Ji n on [0,T C ) D [0,T] we get 

sup ||it(t)|L ^ max( ||u a (t)|L + RJt) ) < +oo . (4.29) 

te[o,T) te[o,T] 

The results ( 13. 8ft and ( 14.29}) are contradictory. 

Step 3. Conclusion of the proof. Eq. ( 14.26}) proved in Step 2 implies 7 = T c ; 
now, the inequality (14.281) of Step 1 coincides with the thesis (I4.27p . □ 

Of course, the control inequalities (14.24)) (14.251) are fulfilled by a function ft n G 
C\[0,T c ), [0,+cx))) such that T c G (0,TJ, and 

= -v^n + (G n V n + K n V n+1 )R n + G n % 2 n + e„ everywhere on [0, T c ), (4.30) 

Bn(0) = 5n • (4.31) 
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4.6 Definition. Eqs. ffJflfy are referred in the sequel as the the control 

equations, or the control Cauchy problem for "R n . 



5 Analytic solutions of the control inequalities 



Throughout this section we consider again the Euler/NS Cauchy problem (13. 5p . 
for a given viscosity v e [0, +oo), external forcing / and initial datum uq as in 
( l3.2p -( pT4"|) . After choosing an approximate solution u & for (I3.5P and determining 
the corresponding estimators, one is faced with the problem of solving the control 
inequalities ( 14. 24ft (I4.25p . or the control equations (I4.30P (I4.3ip . for the unknown 
real function t h-> "Jl n (t). In many applications, such as the one of the next sec- 
tion on the Galerkin approximate solutions, a numerical treatment of the control 
equations is recommended (and the result is generally reliable: (I4.30p (I4.3ip is a 
typically nonstiff Cauchy problem). However, an analytic approach to the con- 
trol equations/inequalities has its own interest, both for theoretical reasons and for 
building simple user-ready criteria. 

In this section we propose an analytic approach for special forms of the ap- 
proximate solution and/or its estimators. For simplicity, throughout the section we 
assume zero external forcing: 



however, many results presented in the section could be extended to the case of 
nonzero /, with suitable assumptions on this function. 

We begin by considering the approximate solution w a := 0; this choice seems to 
be very trivial but, in spite of this, it can be used to obtain nontrivial estimates on 
the time of existence of the exact solution of the Euler/NS Cauchy problem. In the 
NS case v > 0, these estimates include a condition for global existence under a fully 
quantitative norm bound on the initial datum. 

The second case considered in the section is much more general: the approximate 
solution is unspecified, even though a certain form is assumed for its estimators. 

Results from the zero approximate solution. Let us choose for ( 13. 5 p the 

approximate solution 




forcing) 





(5.1) 




w a (t) := for all t . 



(5.2) 
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5.1 Lemma, (i) The differential and datum errors of the zero approximate solu- 
tion are 

e(w a )(t) = for all t G [0, +oo) , w a (0) — u$ = —uo ; (5.3) 
consequently, this approximate solution has the differential and datum estimators 

e n (t) := , 5 n := \\u \\ n , (5.4) 

and the growth estimators 

V n (t) := V n+1 (t) := fortE [0, +oo) . (5.5) 
The control equations l[4-30\ ) l[4-31\ ) with the above estimators take the form 
dH r , 



—— — —v'Jln + G n ^ , 



R n (0) = \\u \\ n . 



(5.6) 



(ii) The Cauchy problem ( 15. 6)) has the solution e R n G C 1 ([0, T c ), [0, +oo)) ; deter- 
mined as follows: 



-oo 



— log 1 



V 



G n \\Uo L 



„ G n ||wo||n 

(intending l/(G n ||wo|| n ) := +°° if u o — 0); 



lfu>0 ! \\u \\n ^ vjG n , 

ifv>0, \\u \\ n > v/G n , 
ifu = 



(5.7) 



3l n (t) :- 



-vt 



1 - G n || Wo 11^(0 
1 - e- ut 



(note that t = lim„_ 



e„(jt) : 
1 - e~ ut 



t 



forte[0,T c 

ifv>0, 
ifv = 



(5.1 



(5.9) 



Proof, (i) Obvious. 

(ii) The Cauchy problem ( 15. 6ft is solved by the quadrature formula 



Uq n 



dr 

G„r 2 — ur 



(5.10) 
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the integral in the left hand side equals (l/z/)log — { n ^J if v > 0, and 

G n - V/\\Uo\\n 

(l/G n )(l/||w ||n — l/%i(t)) if v — 0; some elementary manipulations yield for the 
solution Ji n the expression ( 15.7p -( l5ll . □ 

From the previous lemma and from the main theorem on approximate solutions 
(Proposition 14.5ft . here applied with w a = and 3i n as in 1 15. 8ft . we obtain the 
following result. 

5.2 Proposition. Consider the Cauchy problem A3.5\) for the Euler/NS equations 
with zero external forcing, and any datum u G H™+ 2 ; let u G C([0, T), H" ( | 2 ) D 
C 1 ([0, T), H™ ) be the maximal solution. Define T c and e v as in Eqs. J 5. 7| ) 115. 9\) ; 
then 

T^T C , \\u(t)\\ n ^ S^pT- 7a forte[0,T c ). (5.11) 

1 - G n \\uo\\ n ev{t) 

In particular, due to (5.1), 

T = T C = +oo if \\u \\n < 7^ ; (5.12) 
in this case u is global and, if u > 0, it decays exponentially. 

Other sufficient conditions for global existence. Consider again the Cauchy 
problem ( 13. 5 P with zero external forcing and any datum Uo G H™^ 2 ; let u G 
C([0,T),H£+ 2 ) n C^flp,! 1 ),!^) be the maximal solution. An obvious implication 
of Proposition 15.21 is the following. 

5.3 Corollary. Assume 

K*i)||n<7r for some t x e[0,T) . (5.13) 

Then: 

\\u(ti)\\ e^ 1 ^ 

T = +oo, Kt)lU< - " 1 — — for t e [t u +oo) , (5.14) 

1 - G n \\u{ti)\\ n e„{t - ti) 

with e v as in Eq. A5.9\) . 

Proof. The function u \ [ti,T) is the maximal solution of the Cauchy problem with 
datum u(ti) specified at time ti, rather than at time 0. By Eqs. (15. lip ( I5.12p . with 
an obvious shift in time, we obtain the thesis (I5.14p . □ 

A consequence of the previous result, of more practical use, is the following. 
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5.4 Corollary. Let u & G C([0, T a ), H™+ 2 ) n C^QO, T a ), H™ ) 6e any approximate 
solution of 113. 5\) , with estimators e n , 8 n , V n , T> n+ \; assume that the corresponding 
control inequalities h4-24\ ) M-^ty have a solution Jl n G C 1 ([0, T c ), [0, +oo)) ; with 



T c G (0,T a ]. Finally, assume that 



(V n + tt„)(*i) < for some h G [0,T C ) . (5.15) 



Then: 



Proof. By Proposition 14.51 we have — u a (^)||n ^ for all t G [0, T c ) and, 

in particular, for t = t\\ we further write ||w(ti)||„ ^ ||w a (*i)||n + ll n (^i) ~~ w a (ti)|| n , 
which implies 

\Wt 1 )\\ n < t {V n + 'k n ){t 1 ) . (5.17) 
Eq. ( I5.17P and the assumption (I5.15P gives the inequality 

||«(*l)||n < , 

which has the form (15.13p . By the previous corollary, we have Eq. (I5.14p and this 
result, combined with (I5.17p . gives the thesis f)5.16p . □ 

A general result, under conditions of exponential decay {y > 0) or bound- 
edness (y = 0) for the approximate solution. In this paragraph we exhibit a 
solution of the control equations (I4.30p (I4.3ip . holding for any approximate solution 
whose estimators have certain features. Such features are described in the forth- 
coming Eq. ( I5.18P ; these indicate that the norms of orders n, n + 1 of u & behave 
like e~ ut , while the n-th norm of the differential error behaves like e~ 2ut . (In the NS 
case these are conditions of exponential decay, while in the Euler case they simply 
indicate the boundedness of w a and its error.) 

The assumption that w a behaves like e~ vt corresponds to a rather typical behavior 
of the approximate solutions under zero external forcing; for example, this behavior 
occurs in the case of the Galerkin approximate solutions discussed in the next section 
(see Lemma [6.41 and the subsequent Remark). The differential error of u & typically 
behaves like e~ 2vt when u & is bounded by e~ ut and the differential error depends 
only on the quadratic function 7(u & , u a ); again, this situation occurs in the example 
of the Galerkin approximate solutions with no forcing (see Lemma I6.8p . 

Other approximation methods suggested for the NS equations {y > 0) yield 
approximate solutions with a behavior as in ( I5.18P ; as an example, this happens if 
one assumes no external forcing and takes for u a the truncation to any order of the 
power series solution introduced by Sinai [19]. 
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5.5 Lemma. For any v G [0,+oo), consider for VS. 5\) an approximate solution 
u a G C([0, +00), H™| 2 ) fl C 1 ([0, +00), H™ ). Assume this to have growth and differ- 
ential error estimators of the forms 



V n (t) := D n e~ ut , V n+1 (t) := D n+l e~ ut , £ n (t) := E n e 



-2ut 



(5.18) 



with D n ,D n+ i G (0, +00), E n G [0, +00); furthermore, assume this to have any 
datum error estimator 5 n G [0, +00). From the above constants, let us define 



:{G n D n + K n D n+ i) ; 



(5.19) 



furthermore, assume the "error bound' 



G n E n < F„ 



and define 



W± := F n ± ^Fl - G n E n . 
Finally, consider the function t G [0, +00) 1— >■ e„(i) defined by V5.9\) . and put 

Vn : [0, +00) -> [1, +00) , 7^(t) := ~ ^ M*) . 



(5.20) 
(5.21) 



(5.22) 



Then, (i) (ii) hold. 

(i) The control equations fl^.ffOp jjj4-31\ ) with the above estimators take the form 



~dt 



-v%, + 2F n e- vt H n + G n X + E n e 



-2vt 



^n(O) = 5 n . 



(5.23) 



(ii) The Cauchy problem ( f5. 23\) has the solution H n G C 1 ([0, T c ), [0, +00)), deter- 
mined as follows: 



T c := I 



-00 

--log ( 1 

V 

1 



log 



W+ + GJ n 
W+ - W~ "° W~ + GJn 

W+ + GJ n 



log 



ifu>0, > eW + -w-)/« 

W n -r(j n o n 



i/i/ = ; 



1 W+(W~ + GJ n )rj n (t) - W~(W+ + GJ n ) 



(5.24) 



G n (W+ + GJn) - (W- + GJ^if) 
Proof, (i) Obvious. 

(ii) We write the unknown solution !R n of (I5.23j) as 



e - "* forte [0,T C ). (5.25) 



ft n (t) = Z n {t) e 



-vt 



(5.26) 
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where Z n G C 1 ([0, T c ), [0, +00)) (and T c ) are to be found. Eq. (I5.23P is equivalent 
to the Cauchy problem 

^ = (E n + 2F n Z n + G n Z 2 n ) e~ vt , Z n (0) = 5 n , (5.27) 
at 

with F n as in (15.191) ; this is solved by the quadrature formula 



s n En + 2F n z + G n z 2 jq 



dse~ us . (5.28) 



dse~ us = e v {t) ; (5.29) 
furthermore E n + 2F n z + G n z 2 = (l/G n )(W+ + G n z)(W~ + G n z), which implies 



On the other hand 

ct 







dz 1 



E n + 2F n z + G n z 2 W+ - W n 



4- 1 ri ~~\ Z=Zn{t) 



log 



+ G n z 



(5.30) 

2 = 5,! 



Inserting Eqs. (15. 29ft (I5.30p into Eq. (I5.28P we easily obtain an explicit expression 
for Z n (t) which implies Eq. (I5.25f) for "Jl n {t) = Z n {t)e~ ut . The interval [0,T C ) where 
Z n is well defined is also easily determined, by the same manipulation employed to 
make explicit Z n (t). □ 

From the previous lemma and from the main theorem on approximate solutions 
(Proposition 14. 5j) . we obtain the following result. 

5.6 Proposition. Consider the Euler/NS Cauchy problem Ii3.5\) . and an approxi- 
mate solution w a G C , ([0,+oo),H"+ 2 )nC ,1 ([0,+oo),H™ ) ; assume this to have growth 
and differential error estimators fulfilling all conditions in Lemma 15.51 and define 
T c , Jl n via Eqs. \5.2J^j ( (5. 25\) of the same lemma. Now, consider the maximal 
solution u G C([0,T),H™+ 2 ) n C\[0, T), H£ ) of l[3~5\) ; then 

T^T C , \\u(t)-u a (t)\\^Ji n (t) forte[0,T c ). (5.31) 

(In particular, u is global under the conditions giving T c = +00 in Eq. fi5.24j) -) 

5.7 Remark. One could easily derive a variant of Lemma 15.51 (and Proposition 
15. 6p dealing with the limit case G n E n = F 2 , where W~ = W+. The expressions for 
T c and H n in this case coincide with the ones derivable from Eqs. (I5.24p (I5.25P in 
the limit W~ — > W+. The results in Lemma [57171 (and Proposition 15.21) could also 
be derived from Lemma 15.51 (and Proposition 15.61) . in the limit case E n = and 
D n , D n+1 ->• 0. 
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6 The Galerkin approximate solutions for the Eu- 
ler/NS equations, and their errors 

Throughout this section, we consider a set G with the following features: 

GcZq, G finite , k G G -k G G . (6.1) 

Hereafter we write -< e k >~keG f°r the linear subspace of W made of the sums 
Ylk£G Vkek (which i s > i n fact, contained in C°°). The forthcoming definitions follow 
closely the setting proposed in our previous work [13] for the Galerkin method. 



Galerkin subspaces and projections. We define them as follows. 

6.1 Definition. The Galerkin subspace and projection corresponding to G are, 
respectively: 

W G := D' E0 n -< e k y keG = (^v k e k \v k eC d ,v^ = v. k , k.v k = for all k} ; (6.2) 

keG 

<$ G : D' E0 -> Eg, , v = v ^k := E Vk&k " ( 6 ' 3 ) 





fceZjj k&G 



It is clear that 

Hg, C C°°nD E0 ; A(Hg.) = Eg,; H« c EI™, ^ G (H™) = Hg, for all m G R . (6.4) 
The following result is useful in the sequel. 
6.2 Lemma. Let m,p G R, m ^ p and v G HI E0 . Then, 

\\(l-y G )v\\ m ^ J*- , |G|:= min |*| . (6.5) 



Proof. We have (1 — ( P G )i> = ^2 ke z d \G Vkek '^ thus, 



:i-¥°M|i= E N 2m H 2 = E i=H 2 (6-6) 



A;GZq\G k&Zf\G 



< ( SU P |fc|2p-2m ) E ^l^l^l 2 < 
fcGZ \G I I fc G ZJ*\G 



■INI 2 

£<|2p-2m 11 n P 



□ 



Galerkin approximate solutions. Let us be given z/ G [0, +oo), / G C([0, +oo) 
and Mq G D' E0 . 
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6.3 Definition. The Galerkin approximate solution of the Euler/NS equations 
corresponding to the datum uq, to the forcing f and to the set of modes G is the 
maximal (i.e., nonextendable) solution Uf tU0)G = «g of the following Cauchy problem, 
in the finite dimensional space M G : 



Find uq £ C^O,!©),!!^) such that (6.7) 
^ = is Au G + q3 G P(w G , u G ) + <$ G f , u G (0) = ^ G u . 



Eq. (16. 7j) describes the Cauchy problem for an ODE in the finite-dimensional 
vector space H G ; this relies on the continuous function H G x [0, +oo) — > H G , 
(v, t) i — y uAv + ty G 7{v, v) + *J} G /(t), which is smooth with respect to the variable v. 
So, the standard theory of ODEs in finite dimensional spaces grants existence and 
uniqueness for the solution of (16. 7p . 

The following facts are known (see, e.g., [20J ) ; a proof suitable for the present 
setting is given, for completeness, in Appendix [Bj 

6.4 Lemma, (i) For any v £ [0, +oo) ; / £ C([0, +oo), D' E0 ) and u £ D' E0; the 
maximal solution of problem ^6. 7\ ) is in fact global: 

T G = +oo . (6.8) 

(ii) If the forcing f is identically zero, one has 

\u, mil A = H^oll^ 2 */f = o, t £ [o,+oo), . . 

(Hi) Let v £ [0, +oo). For any f £ C([0, +oo), D E0 ) one has 

\\u G (t)\\ L 2 <C (j|<p G u || L2 + ^ dse vs \\¥ G f(s)\\ L ^e- vt forte[0,+oo) . (6.10) 
In particular, let 

J := / dse vs W G f(s)\\ L2 <+oc; (6.11) 



< (||q3 G %|U 2 + j)e^' /ort £ [0, +oo) . (6.12) 

6.5 Remark. Let v £ [0, +oo). If / = or, more generally, if J < +oo, the 
previous lemma gives a bound of the type ||M G (t)|| L 2 ^ const. e~ vt . This fact, 
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with the equivalence of all norms on the finite dimensional space H^ , implies the 
following: for any real m, 

\\u G {t)\\ m < U m e~ ut for t £ [0, +00) , (6.13) 

with a suitable constant U m £ [0, +00) (depending on the initial datum and on the 
forcing) . 

6.6 Definition. (1) For all k £ Z d , f k £ C([0, +00), C d ) and £ C d are the 

Fourier components of the forcing and the initial datum: f(t) = Ylkez d fhif) e k; 

(ii) For k £ G, 'jck = Ik £ C 1 ([0, +00), C d ) are the Fourier components of ' u G : 

u G (t) = J2^(t)e k . (6.14) 

keG 

(Note the relations fk = f- k , k»f k = 0, and the analogous relations for u ok , Ik-) 
Let us review a well known fact. 

6.7 Proposition. Under the correspondence uq (-> (7^), Eq. ^6. 7\ ) for u G is 

equivalent to the following problem: find a family of functions 7^ £ C 1 ([0, +00), C d ) 
(k £ G) such that 



^ = -v\k\ 2 ^ k - p*y /2 ^2hh»(k - /i)]£ fc 7 fc _ h + f k , 7 fc (0) = u ok (6.15) 



heG 



(intending 7/ c _/ i (£) '■= if k — h (jL G ; recall that is the orthogonal projection of 
C d onto k- 1 ). A family of functions "f k (k £ G) fulfilling Eqs. \6.15\) automatically 
possesses the properties 7^ = 7.^ and k»jk = 0. 

Proof. Clearly, Eq. (16. 7h is equivalent to 

'du G 



dt 



v(Au G ) k + 7(u G , u G )k + fk , (wc)fc(O) = u ok for k £ G ; 



k 



making explicit the above Fourier components via Eq. (I2.23p . we obtain Eqs. (16.151) . 
Let us show that Eqs. (16.151) imply 7^ = 7_ fc and k»^ k = 0. The first implication 
is proved noting that the functions (jk) k e G and (j-^keG are solutions of the same 
Cauchy problem. The second implication follows noting that, for each k £ G, 
(d/ dt){k»^k) = — ^\k\ 2 {k»^k) and (&>7fc)(0) = 0; the solution of this Cauchy problem 
is identically zero. □ 

The Galerkin solutions in the general framework for approximate solu- 
tions. From now on we stick to the framework of the previous sections, i.e.: we 



23 



consider the Sobolev spaces of orders n,n + 1, n + 2, for a fixed n G (d/2 + 1, +oo); 
as in fl3.2|) - fl3.4p . we choose a viscosity v G [0, +oo), a forcing / G C([0, +oo),H" ) 
and an initial datum tt G HI™+ 2 ; we consider the Euler/NS Cauchy problem (13.51) . 
Having fixed a set of modes G as in (16. ip . we regard the corresponding Galerkin 
solution uq as an approximate solution of the Euler/NS Cauchy problem (13. 5p ; our 
aim is to apply the general theory of the previous sections with u SL = uq- 

The desired application requires to give growth and error estimators for uq- Of 
course, we have the tautological growth estimators 



to be used with m = n or m = n + I; these are employed systematically in the 
sequel. Let us pass to the errors of uq and their estimators. 

6.8 Lemma, (i) The Galerkin solution uq has the datum error 




(6.16) 




(6.17) 



keZ d \G 



and its tautological estimator 




(6.18) 



There is a rougher estimator 



u G (0) - u \\ n ^ 5' n 



np • 



Q\p-n 



(6.19) 




(6.20) 



the Fourier representation of the above summands is 




(6.21) 



fcezjj\G 




(6.22) 



k&dG 



24 



dG := (G + G)\(G\J {0}) , p k := -J^yjl ~ h )]^lk-h ■ 

^ n ' heG 

(In the above: G + G := {p + q\ p, q G G}; \ is the usual set-theoretical difference; 
again, we intend ^k-h '■= if k — h G" G.) 
(Hi) The above terms in e(u G ) have norms 



:i-y G )f\\ n = J J2 w 2r w> ( 6 - 23 ) 

fcezg\G 



||(i - <p G mu G , u G ) \\ n = J J2 \ k \ 2n \p k \ 2 ; (6.24) 

these admit the bounds 



k&dG 



ll(l-^)/||n< (6.25) 

||(1 -y G )nu G ,u G )\\ n ^ j^\\u G \\ q \\u G \\ q+1 (6.26) 

where: p is any real number such that p ^ n and f G C([0, +oo), H^ ); q is any real 

number such that q ^ n; K q G (0, +oo) is such that \\y(v,w)\\ q ^ K q \\v 

for all v G H| 0J w G (o/ course ||it G || g = fc S G I^I 2i3 |7a:| 2 j and similarly 

for \\u G \\ q+ i). Thus, a differential error estimator e n of order n for u G is obtained 

setting 



( r.h.s. of ( [Og ) or (gjfflp ) + ( r.h.s. of ( TP3|) or / f031) ) . (6.27) 



Proof, (i) Eqs. (I6.17P f l6T8]) are obvious. Eq. (16. 19|) follows writing ||ttrr(0)— nn||^ = 
||(1 — < P g )m ||„ and using the inequality (I6.5p . 
(ii) We have 

= ^ - " Au G - a>(u G , u G )-f = ¥ G nuG, u G ) + <£ G / - 0>(u G , « G ) - / , 
at 

where the first equality is just the definition of the differential error, and the second 
follows from (j6.7P ; so, Eq. (I6.20p is proved. 

The Fourier representation (I6.2ip for (1 — i $ G )f follows immediately from the 
definition of *}3 G . To derive the representation f)6.22p for (1 — < $ G )'J } (u G , u G ), let 
us consider the Fourier components of u G that we indicate in any case with j k , 
intending lk := for k G Z d \ G. Eq. fl2T23|) gives 



?{u G ,u G ) = 



Pkdk 
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where, for any k, pk is defined following Eq. (I6.22p . On the other hand: for each 
k G Zq, pk is a sum over h G G containing terms of the form 'jk-h', if k ^ G + G, for 
all h G G we have k — h^G (since k — h G G would imply k = (k — h) + h G G + G); 
k — h (jL G implies ^k-h — 0. Due to the above considerations we have pk = for 
k ^ G + G, whence 

7(u G ,u G )= Yl P^k- ( 6 - 28 ) 

ke(G+G)\{0} 

Application of 1 — removes from the above sum the terms with k G G; thus 

(1 - ^ G )T(n G , u G ) = J2 ^ ' ( 6 - 29 ) 

ke(G+G)\(GU{0}) 

in agreement with Eq. (16. 22 p . 

(iii) Eqs. ( I6.23P (I6.24p are straightforward consequences of the Fourier representa- 
tions (JS2IJ (JS23). Eq. (E2SD follows immediately from (153]) . To infer Eq. (15^5]) . 
we note that (16. 5 P and the inequality involving iT g imply 

||(1 -^ G )?(u G ,ua)\\n < j^\muG,u G )\\ q < j^\\u G \\ q \\u G \\ q+1 . 

Finally, Eq. (I6.20p implies 

\\e(u G )\\ ^ \\(l-¥ G )?(u G ,u G )\\ n +\\(l-¥ G )f\\ n ; (6.30) 

binding the two summands in the right hand side via Eqs. (I6.24p or (16. 26 p . ( I6.23P 
or (I6.25P we obtain the estimator e n in (I6.27p . □ 



6.9 Remarks, (i) The exact expression (I6.24p for ||(1 — <! ^ G )7{u Gl u G )\\ n involves 
a sum over the finite set dG; for each k G dG, the term pk in the above sum is itself 
a finite sum, that can be computed explicitly from the Fourier coefficients 7& = 7fc(£) 
of the Galerkin approximate solution (assuming these ones to be known, say, from 
the numerical integration of the system (I6.15P ). However, when G is large the set dG 
is typically very large, and this can make too expensive the computation of the sum 
over dG. In these cases one can use for ||(1 — < Hp G )'J'(u G , u G )\\ n the bound (I6.26P ; this 
requires computation of HmgIIq = V^fceG I^| 29 l7fc| 2 an d of ||wg|| 9 +i, both of them 
less expensive since these sums are over G, rather than dG. 

(ii) We are aware that, if the Galerkin equations (I6.15P are solved numerically by 
some standard method for ODEs, one does not obtain the exact solutions 7fc(£) 
(k G G) but, rather, some approximants whose distance from the 7fc's should be 
estimated on the grounds of the employed numerical scheme. In the application 
presented in the next section, relying on a relatively small set G of modes, the 
intrinsic error in the numerical integration of ( 16 .151) has been regarded as negligible; 
the situation would be different using a much larger set of modes. 
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7 An application of the previous framework for 
the Galerkin method 



A preliminary. In this section we frequently report the results of computations 
performed with MATHEMATICA. A formula like r = a.bcde... must be intended 
as follows: computation of the real number r via MATHEMATICA produces as an 
output a.bcde, followed by other digits not reported for brevity. 

Setting up the problem. Throughout this section, we work in space dimension 

d = 3, (7.1) 

with any viscosity v G [0, +oo). We consider the Euler/NS equations on T 3 with no 
external forcing, in the Sobolev framework of order n = 3. So, the Cauchy problem 
( 13.51) takes the form 

Find u G C([0,T),Wl o ) nC 1 ([0,T),l£ o ) such that (7.2) 

du 

— = uAu + 7(u, u) , u(0) = u . 
dt 

The initial datum uq in H| (in fact, in for any real m) is chosen of this form: 

u = ^2 u okek , (7.3) 

fc=±a,±6,±c 

a := (1,1,0), 6:= (1,0,1), c := (0, 1, 1) ; 
u 0t±a := (2tt) 3 / 2 (1, -1,0), u , ±b := (2tt) 3 / 2 (1, 0, -1) , u ,± c := (2tt) 3 / 2 (0, 1, -1) . 

For v = 0, this initial datum has been considered by Behr, Necas and Wu [2] as 
the origin of a possible blow-up for the Euler equations. More precisely, the authors 
of the cited work try a solution u of the Euler equations in the form of a power 
series u(t) = Y2t^o u i^i where the zero order term corresponds to the initial datum, 
and Ui,u 2 , ... : T 3 — > R 3 are determined recursively. By a means of C++ program, 
the authors compute the terms Ui for i = 1, 35, allowing to construct the partial 
sums Ujv(i) := YliLo u i^ 1 up to = 35. A merely "experimental" analysis of these 
partial sums and of their H| norms brings the authors to conjecture that the exact 
solution u(t) of the Euler Cauchy problem blows up in H| for t — > T*, for some 
T* G (0.32, 0.35). 

For subsequent use, we mention that Eq. f lT.3 j) implies 

||u |U = V37r 3 2 m + 5 (7.4) 
for any real m. This formula gives, for example, 

||uo||i = 77.15..., ||u || 2 = 109.1..., (7.5) 
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||uo|| 3 = 154.3...., \\u \U = 218.2...., ||« || 5 = 308.6... . 

Introducing our approach. In this section we propose a different approach to the 
Euler Cauchy problem of [2], that we apply as well to the NS case; so, we consider 
the problem (17. 2p (17. 3p . for arbitrary v G [0, +oo). 

We refer to the general setting developed in the present paper, using the spaces H£ , 
H™^ 1 , H™+ 2 with 

n = 3 . (7.6) 

This requires, amongst else, the numerical values of two constants K 3 , G 3 fulfilling 
for n = 3 the "basic inequality" (EOl]) and the "Kato inequality" f l2T25|) . Due to the 
computations in [TBJ [TB], these can be taken as follows: 

K 3 = 0.323 , G 3 = 0.438 . (7.7) 

To illustrate our setting, we start with a very elementary result. 

A simple sufficient condition for global existence (and exponential decay). 

According to Proposition 15.2} the solution u of the NS Cauchy problem (17. 2p (17. 3p 
is global and exponentially decaying for t — > +00, if 

^ ^ CslluolU = 67.58... (7.8) 

(in the last passage, we have used the numerical values in (I7.5P (17. 7p for G 3 and 
1 1 ^0 1 1 3) (0)- As shown hereafter, a more refined application of our setting for ap- 
proximate solutions significantly improves the above condition: in fact, in the next 
pages, combining this setting with the Galerkin method we will be able to infer 
global existence for v > 8. 

Going on in our approach: Galerkin approximants. The idea developed in 
the sequel, both in the Euler case [y = 0) and in the NS case {y > 0), is the following: 
to compute numerically the Galerkin approximate solution uq for a suitable set of 
modes G; to construct for it error and growth estimators, in the Sobolev norms of 
orders 3 or 4; to solve numerically the control equations (I4.30p (I4.3ip for an unknown 
function ^3 : [0,T C ) — > [0, +00). After finding the solution of this control problem, 
we can grant on theoretical grounds that the solution u of the Euler/NS Cauchy 
problem (17. 2 p exists at least up to time T c , and that \\u(t) — UG{t)\\ 3 ^ ^{t) f° r 

te[o,T c ). 

2 For completeness, we mention a criterion for NS global existence in H^, (T 3 ), discussed in [14] ; 
this can be written as v ^ ||wo||i/0.407, where uq is an arbitrary initial datum in H^, . With the 
datum uo in (|7.3p (and the value of ||«o|!i in (TZIS])), we conclude that [H] ensures global existence 
in H^ for v ^ 189.5... . By a known regularity theorem about NS equations pUI, this result of 
global existence in H^, also implies global existence in H 3 ; however, the condition v ^ 189.5... 
arising from the H 1 setting of [14] is manifestly weaker than the condition (|7.8p . 
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Our computation has been performed using Mathematica on a PC, with the 
relatively small set of 150 modes 

G:=SU-S; -S := {-k \ k e S} ; (7.9) 
S := {(0, 0, 2), (0, 1, -3), (0, 1, 1), (0, 1, 3), (0, 2, 0), (0, 2, 2), (0, 3, -1), (0, 3, 1), (0, 3, 3), 

(1, -3, -2), (1, -3, 0), (1, -3, 2), (1, -2, -3), (1, -2, -1), (1, -2, 1), (1, -2, 3), 
(1,-1, -2), (1,-1,2), (1,0, -3), (1,0,1), (1,0,3), (1,1, -2), (1,1,0), (1,1,2), (1,2, -3), 
(1, 2,-1), (1, 2, 1), (1, 2, 3), (1, 3, -2), (1, 3, 0), (1, 3, 2), (2, -3, -3), (2, -3, -1), (2, -3, 1), 
(2, -3, 3), (2, -2, -2), (2, -2, 2), (2, -1, -3), (2, -1, -1), (2, -1, 1), (2, -1, 3), (2, 0, 0), 
(2, 0, 2), (2, 1, -3), (2, 1, -1), (2, 1, 1), (2, 1, 3), (2, 2, -2), (2, 2, 0), (2, 3, -3), (2, 3, -1), 
(2, 3, 1), (2, 3, 3), (3, -3, -2), (3, -3, 2), (3, -2, -3), (3, -2, -1), (3, -2, 1), (3, -2, 3), 
(3, -1, -2), (3, -1, 0), (3, -1, 2), (3, 0, -1), (3, 0, 1), (3, 0, 3), (3, 1, -2), (3, 1, 0), (3, 1, 2), 
(3, 2, -3), (3, 2, -1), (3, 2, 1), (3, 2, 3), (3, 3, -2), (3, 3, 0), (3, 3, 2)} . 

The results we present here are somehow provisional; we plan to attack the problem 
by more powerful numerical tools in a future work, using for G a larger set of modes. 

A sketch of the operations to be done. The list of such operations is the 
following: 

(i) First of all, one chooses a value v 6 [0, +oo) for the viscosity, and a finite time 
interval [0, T{) for the numerical computation of the Galerkin solution. 

(ii) The Galerkin solution uq for the set of modes G in ( I7.9P and for the initial datum 
uq is found numerically on the chosen time interval [0, T{). More precisely, one solves 
numerically the system of differential equations (I6.15P for the Fourier components 
(lk)keG of uq, with the initial conditions u k corresponding to Eq. (17. 3p . We recall 
(see Lemma I6.4p that the general theory of the Galerkin solutions for zero external 
forcing ensures the global existence of uq (and its exponential decay, if v > 0); thus, 
from a theoretical viewpoint there is no obstruction to the computation of uq on 
any finite interval [0,T/). 

(iii) We apply to ug on [0, T{) the framework of the present paper for the approximate 
Euler/NS solutions using (we repeat it) the spaces H" , H™^ 1 and EI™+ 2 with n = 3. 

(iv) Some important characters in our approach are the norms 

V n {t) := \\u G (t)\\n = /]T|A;| 2 ™|7 fe (t)| 2 (n = 3,4) , (7.10) 

y fceG 

which can be obtained from the numerical values jk(t) of the Fourier components. 
Further, we need a datum error estimator 63 and a differential error estimator £3. 
On the other hand the datum error is zero in this case, since the initial condition uq 
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is in the subspace M^ spanned by the chosen set (17.91) of modes; thus, we can take 

S 3 = . (7.11) 

As for the differential error estimator, we take the precise expression coming from 
Eqs. (" 16.241) ( 16.27!) (taking into account that the forcing / is zero in this case); these 
yield the expression 

^)=|EW 6 Wi)f (7-12) 

V k£dG 

where dG := (G+G)\(GU{0}) andp fc (t) := -i(2v)- d ' 2 £ fceG [7h(*M*-/0]£*7*-ft(*), 
as in Eq. ( I6.22p . The computation of e 3 (t) following Eq. (I7.12p is rather expen- 
sive: the set dG consists of 929 modes and, for each one of them, one must perform 
the nontrivial computation of Pk{t). (For these reasons, a computation with a set 
of modes much larger than the G in (17. 9p would suggest to replace this e 3 with a 
rougher estimator, coming from Eqs. ( I6.26P ( I6.27P ). 
(v) Now we consider the control equations ( I4.30P ( I4.3ip . taking the form 

dH 

_£ = + (Q 3 V 3 + K 3 V 4 )X 3 + G 3 0q + e 3 , (7.13) 

tt 3 (0)=0, (7.14) 

with K 3 ,G 3 as in Eq.flTTD and V 3 {t), V^t), e 3 (t) as in Eqs. (T7"TUj) (17321 . The 
unknown is a function Jl 3 G C 1 ([0,T c ), [0,+oo)), with < T c ^ Tj; this is deter- 
mined numerically (with a package for self-adaptive integration, allowing to detect 
a possible blow-up of Jl 3 ; in this case, T c is the blow-up time). 

Once we have 3l 3 , the general theory allows to state the following. 

(a) The maximal solution u of the Euler/NS Cauchy problem (17.21) (17.31) has a 
domain containing [0,T C ), and 

\\u{t)-u G {t)\\ 3 ^%{t) for te [0,T C ) (7.15) 

(see Proposition I4.5p . 

(b) Suppose v > 0, and 

(V 3 + ft 3 )(ti) ^ 7^ for some h G [0,T C ) . (7.16) 

G 3 

Then, the solution u of the NS Cauchy problem (I7.2p is global, and 

ll"WII'< i-KS(^- tl ) fOTfe ^ • <" 7 > 

with e u as in Eq. (15. 9p (see Corollary I5.4p . 
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Hereafter, we report the results obtained from the above scheme for some choices of 
v. 

Case v = . First of all, the Galerkin equations ( 16 . 1 5f) have been solved for the set 
of modes (17. 9p on a time interval of length Tj = 2. (This required, approximately, 
15 seconds of CPU time on our machine.) In Figures la, lb and lc we report, as 
examples, the graphs of |7fc(£)| for t G [0,2), in the cases k = (1, 1,0), k = (0,0,2) 
and k = (0,1,-3), respectively. (Of course, we could consider many alternative 
choices, such as plotting the norms |7fc(t)| for other modes, or the real parts of the 
components 7£(i) (r = 1, 2, 3), or the imaginary parts of the same components.) As 
a supplementary information, we mention that the graphs of |7fc(£)| are identical in 
the three cases k = (1, 1, 0), (1, 0, 1), (0, 1, 1)). 

After computing numerically all the 7fc(t) (for k G G), one obtains from Eqs. 
(17.1 Op and ( I7.12p the norms V n [t) := ||w G (i)|| n (n = 3,4) and the differential error 
estimator es(t), for t G [0,2); Figures Id and le report the graphs of T)$(t) and 
63(f) (to give some more detail, we mention that the functions T> 3 (t), V±(t) and e 3 {t) 
have been computed for a set of sample values of t G [0,2), and then interpolated 
by means of the MATHEMATICA algorithms. The computation of e 3 (t) at each 
sample value t is rather expensive, for it involves a sum on the large set dG; the 
CPU time is about 15 seconds for each t and, for this reason, we have used only 30 
sample values). 

The final step is the numerical solution of the Cauchy problem (I4.30P (I4.3ip 
for the unknown function ft 3 (t) (a task performed almost instantaneously by our 
PC). The MATHEMATICA self-adaptive routines for ODEs indicate a divergence 
of R 3 (t) for t -»■ T c , with T c = 0.06666... . Figure If gives the graph of CR 3 (t) for 
t G [0, T c ). The conclusion of these computations is the following: the solution u 
of the Euler Cauchy problem (17. 2p (17.31) has a domain containing [0, T c ), and its 
distance from uq is bounded by Ji 3 on this interval. 

Some remarks on the i> = case. We have already mentioned that the blow-up 
of u(t) is conjectured in [2] for t — » G (0.32, 0.35). The basis of this conjecture is 
an experimental analysis of the first 35 terms in the power series (in time) solving 
formally the Cauchy problem; in principle, this analysis does not even prove exis- 
tence of u(t) at any specified time t < T*. Our lower bound T c = 0.06666... for the 
interval of existence of u is about 1/5 of the blow-up time T* suggested by [2], but it 
relies on an analytic theory for approximate solutions, their errors, etc., summarized 
by Proposition I4.5( in this sense, it is theoretically grounded. 

Perhaps, our method could give a sensibly larger lower bound T c on the time of 
existence, when implemented with a set G of Galerkin modes much larger than (17.91) . 
Alternatively, one could apply our theoretical framework using as an approximate 
solution the partial sum un 

(*) : = Ef=i^\ with N = 35 as in [2] or with a larger 
N. Both tasks require much more expensive numerical computations (to be done 
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with more appropriate hardware and software): we plan to do this elsewhere. One 
cannot exclude that an attack with more powerful devices could finally result into 
a theoretically grounded lower bound T c larger than the suspected blow-up time T* 
of [2]; however, at present this is just a hope. 

Cases v = 3 and v = 7. We use again the Galerkin solution uq with G as in 
f)7.9p . Due to the positivity of u, all components Jk(t) °f the Galerkin solution 
decay exponentially for t — > +oo (recall Lemma I6.4p ; the same happens of the 
norms T> 3 (t), T>±(t) and of the differential error estimator e 3 (t), which are essential 
objects for our purposes. 

The system (I6.15P for the Galerkin components 7fc(t) has been solved numerically 
on a time interval of length Tj = 1 (which required a CPU time of about 15 seconds 
for v — 3, and 25 seconds for v = 7). Subsequently, V 3 (t), T>±(t) and e 3 {t) have 
been computed from the components 7fc(t) and Eqs. (I7.10p (I7.12p (indeed, some 
interpolation has been done as in the case v = 0; as in that case, for the computation 
of 63(f) we have used only 30 sample values of t in [0, 1), with a CPU time of about 
15 seconds for each one). 

For both the above values of u, the final step has been the (very fast) numer- 
ical solution of the Cauchy problem (I4.30p ( 14 . 3 1|) for the unknown function "Jl 3 (t). 
According to the MATHEMATICA routines for ODEs, Jl 3 (t) diverges for t T c , 
where T c = 0.09025... for v = 3, and T c = 0.2386... for v = 7. We repeat that these 
results grant existence on a domain D [0,T C ) for the solution u of the NS Cauchy 
problem (17. 2p (17. 3p . and the bound (I7.15P on this interval. 

As examples, in Figures 2a-2f we have reported some details on computations for 
v = 7. More precisely, Figures 2a, 2b and 2c are the graphs of |7fc(£)| for t 6 [0, 1), 
in the cases k = (1,1,0), k = (0,0,2) and k = (0,1,-3), respectively. (In fact, 
the graphs of |7fc(t)| for k = (1,0,1) and k = (0,1,1) are identical to the graph 
of the case k = (1,1,0).) Figures 2d and 2e are the graphs of T> 3 (t) and e 3 (i), 
for t G [0,1). Figure 2f gives the graph of ^(t) that, as anticipated, diverges for 
t^T c = 0.2386... . 

Case v = 8. Again, we have used the Galerkin solution uq with G as in (17. 9p . All 
the components 7fe(t), the norms T> 3 (t), T>^(t) and the differential error estimator 
e 3 (t) decay exponentially for t — > +00. The system (I6.15P for the Galerkin compo- 
nents 7fc(t) has been solved numerically on a time interval of length Tj = 1 (which 
required a CPU time of about 25 seconds). The forthcoming Figures 3a, 3b, 3c are 
the graphs of |7 fc (t)| for t £ [0,1), k = (1,1,0), k = (0,0,2) and k = (0,1,-3), 
respectively. Subsequently, the norms V n (t) [n = 3,4) and the error e 3 (t) have been 
computed from the components 7^(t) and Eqs. (I7.10p (I7.12p (making some inter- 
polations, as in the previous cases). Figures 2d and 2e are the graphs of T> 3 (t) and 
e 3 (t), for t £ [0,1). 

The final step has been the (very fast) numerical solution of the Cauchy problem 
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f)4.30p (I4.3ip for the unknown function 3? 3 (t). Differently from all the previous 
cases, the numerical solution 3l 3 (i) determined by MATHEMATICA is defined on 
the whole interval [0, 1); its graph is reported in Figure 3f which suggests, via some 
extrapolation, that 3? 3 (£) should be defined for all t G [0, +oo), with $ 3 (i) — > + for 
t — > +oo; of course, this would imply global existence for the solution u of the NS 
Cauchy problem (Q (Ojl . 

However, global existence of u can even be inferred without extrapolating the 
behavior of 3l 3 outside the interval [0, 1). In fact, global existence is granted if the 
condition (17.161) (V 3 + ^3) (ti) ^ vjG 3 holds at any instant t\ > 0; in the present 
case vjG 3 = 18.26..., and the numerical computation performed in the interval [0, 1) 
shows that (I7.16P is satisfied for any ti G [0.1567..., 1). In conclusion, we can take for 
granted that we have global existence for the solution u of the NS Cauchy problem. 
Of course, \\u(t) — w a (£)||3 is bounded by the numerically computed function "Jl 3 (t), 
for t G [0,1). After choosing a t\ G [0.1567, +00), we obtain as well a bound of 
the form (I5.14p for ||-u(t)|| 3 , which also implies ||w(t)||3 to vanish exponentially for 
t -> +00. 

For example, let us choose £1 = 0.9. We have D 3 (0.9) = 8.580... x 10~ 5 , ft 3 (0.9) = 
0.06100...; with the already known value G 3 = 0.438, Eq. (I5.14p gives 

0fil4p- 8 (*-°- 9 ) 

Mt) h < 1 + 0.00335 e -8(t-o-9) < °- 0614 e " ( "°- 9) for * G ^ +°°) " ( 7 ' 18 ) 

(Here the first inequality follows directly from (I5.14p . recalling that e§(t) = (1 — 
e _8t )/8; the second inequality is obvious.) 

Summary of the previous results, and final comments. Our method to treat 
the Galerkin approximant ug, with G as in (17. 9p . has given the following results 
about the Euler/NS Cauchy problem (17.21) (17. 3p . 

a) v = 0, 3, 7: we can grant existence of the solution u of (17. 2p (17. 3p on an interval 
containing [0,T C ), with T c = 0.06666,0.09025,0.2386, respectively, for these three 
choices of v. We have ||w(£) — u G (t)\\ 3 ^ CR 3 (t) for t G [0,T C ), where H 3 is computed 
numerically solving (I4.30p (I4.3ip (for v = and v = 7, the graph of 3£ 3 is reported 
in Figures If and 2f). 

b) v — 8: we can grant global existence for the solution u of (I7.2p (I7.3p . For 
t G [0, 1) we have a bound \\u{t) — mg(^)||3 ^ 3? 3 (t), with "Ji 3 obtained again from the 
numerical solution of ( 14.301) (14.3 ip ; the graph of this function is reported in Figure 
3f. For t G [0.9, +00) we have a bound of the form (I7.18P on ||it(£)|| 3 , decaying 
exponentially for large t. 

By extrapolation, we are led to conjecture that results similar to (a) should be 
obtained for all v G [0,i> cr ), while results similar to (b) should be obtained for all 
v G [u cr ,+oo), for some v CT G (7,8). (From a qualitative viewpoint, this is just the 
behavior described by Lemma Iq"T51 on the analytical solution of the control inequality. 
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However, here we are applying the control equalities with the tautological growth and 
error estimators, given directly by the numerical solution of the Galerkin equations; 
these are more precise than the simple analytical estimators considered in the cited 
lemma.) 
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Figure la. v = 0. Graph of |7fc(t)| when k 
(1,1,0), forte [0,2). 




Figure lb. v = 0. Graph of |7fc(t)| when k 
(0,0,2), forte [0,2). 




1.5 



Figure lc. v = 0. Graph of |7fc(i)| when A; 
(0,1,-3), for t e [0,2). 
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Figure Id. v = 0. Graph of D 3 (i) = 
||u G (t)|| 3 , for t e [0,2). One has D 3 (0) = 
154.3..., 2? 3 (0.02) = 156.4..., X> 3 (0.04) = 162.7..., 
P 3 (0.06) = 172.7... . 




Figure le. v = 0. Graph of e 3 (t), for f e [0, 2). 
One has e 3 (0.02) = 9.027..., e 3 (0.04) = 36.17..., 
e 3 (0.06) = 81.65... . 
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Figure If. v = 0. Graph of 3? 3 (t); this func- 
tion diverges for t — > T c = 0.06666.... One has 
ft 3 (0) = 0, 3? 3 (0.02) = 0.1439..., 3? 3 (0.04) = 
4.685..., 3? 3 (0.06) = 182.3... . 



35 




0.2 



0.4 



0.6 



0.8 



Figure 2a. v = 7. Graph of |7fc(t)| when k 
(1,1,0), forte [0,1). 
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Figure 2b. v = 7. Graph of |7fc(t)| when k 
(0,0,2), forte [0,1). 
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Figure 2c. v = 7. Graph of |7fc(t)| when k 
(0,1,-3), for t e [0,1). 




0.4 0.6 0.8 1 

Figure 2d. v = 7. Graph of £> 3 (t), for 
t e [0, 1). One has X> 3 (0) = 154.3..., X> 3 (0.025) = 
109.5..., r> 3 (0.07) = 58.40..., X> 3 (0.15) = 18.89..., 
X> 3 (0.23) = 6.153... . 
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0.2 0.4 0.6 0.8 1 

Figure 2e. v = 7. Graph of e 3 (t), for 
t e [0,1). One has e 3 (0) = 0, e 3 (0.025) = 
2.024..., e 3 (0.07) = 0.5788..., e 3 (0.15) = 0.01089..., 
e 3 (0.23) = 1.338... x 10" 4 . 
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Figure 2f. v = 7. Graph of ft 3 (t); this 
function diverges for t — >• T c = 0.2386.... One 
has % 3 (0) = 0, ^ 3 (0.07) = 2.096..., #3(0.15) = 
20.90..., # 3 (0.23) = 265.0... . 
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Figure 3a. v = 8. Graph of |7fe(t)| when k 
(1,1,0), forte [0,1). 
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Figure 3b. v = 8. Graph of |7fc(i)| when k 
(0,0,2), forte [0,1). 
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Figure 3c. v = 8. Graph of |7fc(i)| when fc 
(0,1,-3), for t e [0,1). 
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Figure 3d. v = 8. Graph of D 3 (t), for 
< e [0,1). One has P 3 (0) = 154.3..., X> 3 (0.2) = 
6.280..., X> 3 (0.4) = 0.2559..., X> 3 (0.6) = 0.01043..., 
X> 3 (0.9) = 8.580... x 10~ 5 ... . 
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Figure 3e. v = 8. Graph of e 3 (t), for t 
[0,1). One has e 3 (0) = 0, e 3 (0.023) = 1.553. 
e 3 (0.07) = 0.2844..., e 3 (0.15) = 0.002638. 
e 3 (0.23) = 1.593... x lO" 5 . 




.2 0.4 0.6 

Figure 3f. v = 8. Graph of ft 3 (t), for t € [0, 1). 
One has IR 3 (0) = 0, 3? 3 (0.2) = 6.435..., £ 3 (0.4) = 
2.787... # 3 (0.6) = 0.6503..., 3? 3 (0.9) = 0.06100.... 
The condition (|7. 16|) for global existence of the 
NS Cauchy problem, i.e., (V 3 +$. 3 ){t 1 ) v/G 3 , 
is satisfied for any t\ e [0.1567..., 1) . 
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A Some comparison lemmas of the Caplygin type 

Let 7 G (0, +00]; suppose we are given 
ip G C(R x [0,T),R), (s,t) H> ip(s,t) such that ^GC(Rx[0,T),R); (A.l) 

s G R . (A.2) 
Under the above assumptions, the Cauchy problem 

^(t)=V(S(*),t), 8(0) = « . (A.3) 

has a unique maximal (i.e., nonextendable) solution S G C 1 ([0, T§), R). (Of course 
Tg G (0, T]; later on we give conditions under which T§ = 7.) 

The following is a known result, of the Caplygin type: see [9] or [TT] . 



A.l Lemma. Suppose there is a function W G C([0,T),R) stzc/i t/jat 



f/t ■(*) /orte[0,0% W(0)<s o (A.4) 

fwift d+/dt the right lower Dini derivative, see Eq. A2.2\) ). Then 

W(t) <S(t) /ortG [0,T S ) . (A.5) 
A straightforward variant of the previous result is the following. 
A.2 Lemma. Suppose there is a function 31 G C([0,T),R) suc/i that 

— (t) ^ VW),*) [O,^), ^(0) ^ s , (A.6) 

(with d + /dt the right upper Dini derivative, see Eq. A2.3\) ). Then 

&(*)>§(*) forte[0,T s ) . (A.7) 

Proof. We consider the function —Jl, recalling that d+(—3V)/dt = —d + H/dt; from 
( 1A.6I) we infer 



d+ ( ^ (t) < -V»(3l(t),t) for t G [0,T), -01(0) < -s • (A.8) 

On the other hand, the function — S G C 1 ([0, T§), R) is the maximal solution of the 
Cauchy problem 

t^(i) = -^(S(t), t) , -S(0) = -so . (A.9) 
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Therefore, Lemma IA.1I with W replaced by — 31 (and other obvious substitutions) 
gives 

^ -S(t) for t G [0,T S ) , 
yielding the thesis ( IA.7j) . □ 



A. 3 Lemma. Suppose there are functions W, 31 G C([0,T),R) fulfilling Eqs. 
IjA.J^ , liA.6]) respectively. Then 

T s = T , W(t) ^ S(t) s$ K(t) /or t G [0, T) . (A.10) 

Proof. From Lemmas IA. II and IA.2I we infer 

W(t) ^ §(t) ^ R(t) for t G [0,T§) ; (A. 11) 

we claim that this inequality implies 

T S = T. (A.12) 

In fact, whenever T§ is finite, the nonextendability of S beyond this time implies 
that § is unbounded in any left neighborhood of T s ~; on the other hand, if it were 
T s < T, Eq. flA.llj) would ensure the boundedness of § on [0, T§). 



Now, the combination of (IA.11I) and (IA.12I) gives the thesis flA.lOj) . □ 
If we forget the solution S of the Cauchy problem flA.3j) . the inequality flA.lOj) 



becomes W(t) K(t) for t G [0,T); this is just the statement in Eq. fl423|) of 
Lemma \4.'S\ which is thus recognized as a weaker version of Lemma IA.31 



B Appendix. Proof of Lemma [674] on the Galerkin 
appr oximant s 

Let v G [0,+oo), / G C([0, +00), D' E0 ) and u G D' E0 . In the finite dimensional 
space H E0 , we consider the maximal solution uq of the Cauchy problem (16.71) . here 
reproduced for the reader's convenience: 

Find u G G C\[0,T G ),M.g) such that 

^ = is Au G + ^ G nu G , u G ) + f , u G (0) = ^ G u . 
at 

For the moment, the maximal time of existence T G is not known; one of our aims 
is to show that T G = +00. Let us proceed to prove this claim and all the other 
statements of Lemma 16. 4( our arguments are divided in several steps. 
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Step 1 ("energy balance equation"). Everywhere on [0,T G ), one has 

j t |M|| 2 = 2v (Au G \u G ) L 2 + 2(¥ G f\u G ) L 2 . (B.l) 

In fact, 

sMi.^M* < B - 2 > 

= 2z/(A Mg | Mg ) l2 + 2(¥ G ?(u G , u G )\u G ) L 2 + 2(^ G /| MG ) L2 . 
On the other hand, 

(¥ G ?(u G ,u G )\u G ) L 2 = (7(u G ,u G )\u G ) L 2 = ; (B.3) 

in the first passage above, the projection on M G Q is omitted since we are taking 
the inner product with u G = u G (t) G HI G ; the second passage reflects the well known 
identity (7(v,w)\w) ^2 = 0, holding for all zero mean, divergence free vector fields 
v, w on T d sufficiently regular to give meaning to the indicated inner product (see, 
e.g., [IS], Lemma 2.3). Eqs. (lR2j) flB~3|) yield the thesis dB~Tj) . 
Step 2. Everywhere on [0,T G ), one has 

j t \\u G \\ 2 L 2 < -2z/|| ? i G ||i 2 +2||^ G /|| i2 || MG || L2 . (B.4) 

This follows from the energy equation ( IB.lj) . combined with the inequalities 

(Au G \u G } L 2 < -\\u G \\l 2 , (¥ G f\u G } L 2 ^ ||^ G /|| L2 ||n G || L2 . (B.5) 

The first of these relations follows using ( I2.19p . with n = 0; the second one is an 
application of the Cauchy-Schwarz inequality for ( | ) L 2. 
Step 3. Everywhere on [0,T G ), one has 

^ 1Mb ^ -HMb + W G f\\ L * , (B.6) 

with d + /dt the upper Dini derivative. Let us first prove Eq. (1B.6|) at a time to £ 
[0,T G ) such that u G (t ) ^ 0. In this case, ||it G (t)||,L 2 7^ for all t in an interval / 
containing to] in this interval the function t \-¥ ||m g ||l 2 is C 1 , and we have 

4: \WgWl* = ^: hah? = 0| , \ ||«G||ia < -^Ikoll^ + ||^ G /|U2 , 

dt dt 2\\u G \\ L 2 dt 

the last passage following from ( 1B.4I) . Now, consider an instant to such that u G (to) = 
0. Due to a general result already mentioned (see the comments before Eq. ( 14. 18ft ). 
we can write 

d+ II II ^ \\ duG (+ Ml 

dt t IMIl2 ^ n^* 2 ; 
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on the other hand, Eq. (16.71) and the assumption ug(^o) = give (duc/dt)(t ) 
y G f(to),so 

IMIl 2 < W G f(t Q )\\ L2 ; 

to 



dt 



this is just the thesis (1B.6|) . at the instant under consideration. 
Step 4. One has 

\\u G (t)\\ L 2 <: ^ G u4 v > + j\se us \\^ G f(s)\\ L ,)e- vt for t G [0,T G ) . (B.7) 

Let us consider the continuous function t G [0, Tq) ^ \\uG{t)\\ip] due to Step 3 and 
to Ug(0) = ^P G wo, we have 

^\\u G (t)\\L* < il>(\\u G (t)\\ L 2,t) , \\u G (t)\\ L 2 = \\¥ G u \\ L 2 (B.8) 

where 

V> : R x [0, +oo) ->■ R , VM) : = -vs + \\ < # G f(t)\\ L 2 . (B.9) 
On the other hand, the Cauchy problem 

^(t)=V(S(t),t) , 8(0) = H^olUa (B.10) 
has the global solution 

S(t) := (\\¥ G u4 L 2 + ds e» s \\¥ G f{s)\\ L 2^e- vt for t G [0, +oo) . (B.ll) 

So, by the comparison results reviewed in the previous Appendix (see, in particular, 
Lemma [A. 2 1) , we have 



\u G {t)\\ L ^Ht) fortG[0,T« 



G 



this is just the thesis (IB. 71) . 

Step 5. One has Tq = +oo (which justifies statement ( 16.8)) in the lemma under 
proof). If it were Tq < +oo, we would have lim sup^j,- ||wg(^)||l 2 — +oo; this 
would be contradicted by Eq. (IB. 7ft . implying sup t£ r 0jTo ) ||ug(0IU 2 ^ II < P G ' u o||l 2 + 
/ TG ^e-||^ G /( S )|| L2 <+oo. 

Step 6. Proof of statement A6.9\) in the lemma, under the assumption of zero forcing. 
We use the results of Steps 1 and 4, putting therein Tq = +oo. For / = and v — 0, 
the energy balance equation (IB.lj) takes the form [d/dt) ||"Ug||^2 — 0; thus, 



\u G {t)\\ L 2 = \\U G {Q)\\ L 2 = H^MolU 2 
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for all t E [0, +00). For / = and v > 0, Eq. flRTj) gives 



\\u G (t)\\ L , < ||^ G n || L2 e-^ , 

again for t G [0, +00). The above two equations correspond to the content of (16.91) . 
Step 7. Proof of statements ( 16'. l(J\) - ftKT2\) in the lemma. Statement (16.101) is just Eq. 
(IB. 7j) . with T G = +00. Having proved (16.101) . statements (16.111) (16.121) are obvious. 
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